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1 Introduction 

In [6] we initiated a study of conformal mappings to a gearlike domain with 
a single tooth: a starlike open set in the complex plane bounded by arcs of 
two circles centered at the origin and segments of two lines passing through 
the origin. The approach was to examine the Schwarzian derivative of such 
a mapping (a rational function of a particular form). Relationships were 
found between the auxiliary parameters in the Schwarzian derivative and 
the geometry of the gearlike domain, and the set of parameters producing a 
univalent mapping (“region of gearlikeness”) was determined. 

Using the theoretical results in [6], we focus here on the computational 
aspects of the conformal mappings to one-tooth gear domains from both a 
disk and a rectangle. We begin in section [2] with the relationship between 
the parameters in the Schwarzian derivatives for these two cases. In section [3] 
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we work out the transformations or normalizations necessary for the solution 
of the Schwarzian equation to produce a gear (rather than just a “pregear” 
in the terminology of [6]). In section 0] we describe a variety of approaches 
for carrying out the computations, and report some numerical results in 
section |5l This includes a “map” of the internal structure of the region of 
gear likeness. 

Finally, in Section [7] we give three applications of our results on gear 
mappings: evaluation of a singular integral, mapping to the complement of 
an annular rectangle, and symmetric multitooth domains. 

An Appendix contains some background material on one of the numerical 
methods employed. 


2 One-tooth gear domains 

General facts stated below concerning conformal mapping, Schwarzian deriva¬ 
tives, etc., may be found in many standard texts such as [3 031 m [20]. A 
one-tooth gear domain is a topological quadrilateral in the complex plane 
which is the union of an open disk and a sector of a concentric disk of larger 
radius. Thus up to an affine equivalence az + h, a. one-tooth gear domain is 
the standard domain with vertices at wi = W 2 = e*'^, wz = 

W 4 = where the edges [tci, W 2 ] and [tca, are straight segments called 

the tooth edges of the gear, and the edges from W 2 to Wz and from to Wi 
are arcs of the circles {|w| = 1 } and {|tn| = (3} subtending angles in the 
ranges 7 < arg w <271 — ^ and —7 < arg w < 7 respectively. We say that (3 
is the gear ratio and 7 is the gear angle. 

From our perspective, the hrst thing one needs is to have available is 
the Schwarzian derivative Sf = (/"//T “ /f'Y of fh® conformal 

mapping. When Sf is known, a unique solution / is determined by the triple 
J/( 0 ), where the 2-jet of / at any point z is dehned to be 

= ( 1 ) 
as long as we assure that f'{0) 7 ^ 0 . 

As we note in Proposition 2.1 of [5], if /(O) is the tooth center (the 
common center of the concentric circles containing the circular boundary 


2 


arcs of the gear) and /'(O) = 1, we have 

J/(0) = (0,1, 2 (cost 2 — costi)) (2) 

where 0 < < ^2 < and and are the prevertices. However, 

when we solve for / given Sf and with the normalization ([2]), /(O) does not 
generally turn out to be the tooth center. Unless we know precisely what 2- 
jet to use, the solution / will be a Mobius transformation of a one-tooth gear 
domain, which we term a pregear (see Figure |2]). Although pregears are less 
rigid objects than gears, they have some very restrictive properties. Note that 
the tooth edges of any pregear may be uniquely identihed by having different 
interior angles at their endpoints. Applying Euclidean transformations, we 
restrict the discussion to pregears which are symmetric in M and have no 
vertices on M. In [S] the following was shown. 

Proposition 2.1 (a) Let D he a circular quadrilateral with the above sym¬ 
metries, having two interior angles equal to 7i/2 and two interior angles equal 
to 37r/2. Assume that one tooth edge of D lies in the upper and the other in 
the lower half-plane. Then D is a pregear if and only if the full circles , 
C~ containing the tooth edges intersect in two points. 

(b) Let D be a pregear. Then D is a gear if and only if its tooth edges are 
straight, or equivalently, if the non-tooth edges are arcs of concentric circles. 


2.1 Conformal mapping from disk to a one-tooth gear 


Let D = {z G C: \z\ < 1} denote the unit disk, and let /:D —)■ Gi 3 ^.y be 
a conformal mapping. Suppose the prevertices Zi = f~^{wi) are located at 
points of the form and for 0 < U < ^2 < tt. The expression 

for the Schwarzian derivative iS/- of / as a rational function Sf = Rti,t 2 ,\ 
terms of ti, ti, and an auxiliary parameter A was worked out explicitly in [ 6 ]. 
There it was also noted that by means of precomposition of / with a Mobius 
transformation 


Tfz) 


z — q 
1 — qz 


( 3 ) 


(—1 < q < 1) which leaves D invariant, the prevertices can be symmetrized 
with respect to the imaginary axis. Then (writing now / in place of / o Tg) 
the Schwarzian derivative takes the form 


Sf = Rt,x 


( 4 ) 
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where 


1 

2 





(5) 


and 


V’o,t(2:) 


(sin^ t) {z‘^ — (16 cos t)z^ + {4 + 2 cos 2t)z‘^ — (16 cos + 1) 
2 ( 2 :^ — (2cos2t)2:2 + 1)2 


AAA 


—8 cos t 


( 6 ) 


z"^ — (2 cos2t)z‘^ + 1 


One advantage of the symmetrized form is that the conformal M{t) of the 


gear domain is easily related to the single parameter t via an elliptic integral, 
as we will see in a moment. In [ 6 l section 3] this was exploited to derive 
qnalitative relationships between the pair of geometric parameters /9 ,7 that 
prescribe the gear and the mapping parameters t, A. In this paper we are 
more concerned with qnantitative information concerning this relationship. 

2.2 Conformal mapping from rectangle to gear 

In the stndy of circnlar qnadrilaterals with two symmetries in [5], mappings 
from a rectangle to the quadrilateral were investigated as well as mappings 
from a disk. For numerical work such mappings have certain advantages 
over mapping from the disk, among them the fact that certain Schwarzian 
derivatives are real on the boundary (or on horizontal or vertical sections of 
the rectangle). The relationships between conformal mappings to a gear from 
a disk and from a rectangle given below are somewhat more complicated than 
in the situation discussed in [5] but follow the same general line of reasoning. 

As in [S], we let E denote the elliptic integral 



dz 


(7) 


0 AA~A){z + e AA 4- eAA — e **) ’ 


This is a particular case of a Schwarz-Christoffel mapping [7] , and the image 
R = £^(D) = [—a;i,a;i] x [—Ima; 2 , Ima; 2 ] is a rectangle centered at the 
origin. Here cui, U 2 are half-periods of the related Weierstass p-function, 
with oji > 0, \muJ 2 = |ci; 2 |, and r = UJ 2 / 0 J 1 . Thus the composition / o E~^ 
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is a vertex-preserving conformal mapping from a rectangle to the gear G/ 3 , 7 , 
whose conformal module [TS] is by dehnition the ratio 


, , Im Eli) 

of the sides of the rectangle. 

As customary, write e* = p^cji), i = 1,2,3 where = uji + U2- These 
values are real and satisfy 62 < 63 < Ci. Also as in [S] for dehniteness we 
will have cji, 002 normalized so that ei — 62 = 4, so the periods are in fact 
uniquely determined by their ratio r. The change of variable 


maps D onto the subrectangle 

Ro = {C- -y < ReC < y, -Imy < ImC < Imy} (10) 

of i?, and are mapped to the vertices of Rq. We dehne 

= -4 (p(c + + P(C + (11) 

where fi is real (Figure [I]). It is straightforward to verify, using standard 
properties of Weierstrass elliptic functions |25] that 9 ?T,/i(C) R real when C ^ 
and also when ( G ORq (although not for imaginary ( in general). The 
composition 

^(C) = /of;-'(2C) (12) 

is our mapping g: Rq —?• G/ 3 , 7 . 


Proposition 2.2 Let G = G/ 3,7 he a gear domain and let g\ Rq ^ G he a 
conformal mapping which respects vertices, taking the horizontal edges of Rq 
to the tooth edges of G and the left vertical edge of Rq to the inner circum¬ 
ference of G. Then the Schwarzian derivative of g is egual to Sg = (pr^g for 
some values of r and p. 

Proof. The argument is quite similar to that of O eq. (12)], so we will only 
mention the salient points. One extends g by reflection across the edges of 
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Rq and notes that the composition of reflections along opposite edges in the 
C-plane is a translation of magnitude 2ui or 2 u 2 - By the Chain Rule for 
Schwarzian derivatives, one sees that Sg{() is a function which has periods 
of both of these magnitudes. Due to the fact that g sends the right angles 
at the vertices 003/2 = (cui + oj 2)/2 and 0 J 3 to the right interior angles of G, 
a calculation of the series expansion centered at these vertices shows that g 
is holomorphic at these vertices. In contrast, due to the angles of 37r/2 at 
the remaining vertices of G, the Schwarzian derivative Sg has double poles 
at the left vertices of Rq, with singularities 

-4 -4 

(C+W^’ (C + ^/2)2' 

One checks that the function given explicitly by flTTll is an elliptic func¬ 
tion with the same singularities, so the difference is a constant, which is then 
seen to be real due to the symmetry. □ 

The parameter /i plays a role similar to that of the parameter of the 
same name in [5], just as does A for our mappings of the disk, although the 
Schwarzians are now different functions in the present work. 



Figure 1: \(pT,ti\ for Schwarzian derivative flTTl) 
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2.3 Relation between disk and rectangle Schwarzians 


Proposition 2.3 Let 0 < t < 7 r/ 2 , let r/i > 0, and let /: D —>■ C and 
g: Rq ^ C be conformal mappings with Schwarzian derivatives 

Sf{z) = = L>rA0 

respectively. Suppose that g{C) = f{z) where ^ holds. Then r/i = M{t) 
and 

3 + cos2t 

/i=16AcostH -. (13) 

6 

Proof. The relation between t and r is evident by the definition of conformal 
module. Define 

Tl{z) = p{C+ -^-), (p 2 {z) = p{C+ -^-), 

SO 9?r,o(C) = ~4((^i(z) + (p 2 {z)). From the symmetries and the well-known 
fact [25] that p: (C mod {2n;i, 2u2}) —)• CU{oo} is a 2-to-l covering branched 
at the vertices of R, we see that pi and ip2 map D onto the upper and lower 
half planes, respectively. By the Reflection Principle, it can be deduced that 
(pi and (p 2 are Mobius transformations. Further, they send the quadruple 
(e**, —e“**) to 

(ea, 02, oo, Cl), (02, 03, Cl, oo) 


respectively. Note that cross ratios of the above quartets are equal due to 
the relation 01 - 1-02 + 63 = 0. From this it follows that ip 2 = ( 62+1 + 0 ^ + 
0 i 03 )/(+i — 62 ), and then a straightforward calculation verifies that 


Z — 

ipi{z) = 03 +z(0i -03)(cott)-, 

z + 

_ jf 

Z 

T 2 {z) = 03 -i( 0 i - 03 )(cott) -^- - 

Z + 0 “** 

We now substitute these formulas in flTT]) to obtain a formula for +t,h(C) as 
a function of 2 ;, 




(| cos 2t)z^ + (y cos 2t cos t — 8 cos + | cos 2t 
z^ + (2 cos 2 ^) 2 ; + 1 


+ 4p,. 


(14) 
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The Chain Rule for Schwarzian derivatives says that 



and the terms involving E are easily evaluated from the definition ([7]) as in 
[5]. After simplification, we find that 


cqd + ciz + C2Z^ + + CiZ‘^ 

3A2 



Sf{z) 


where cq = C 4 = —2cos2f, ci = C 3 = 12(cos3t — cost), C 2 = 5 — cos4f, 
A = z'^ — (2cos2t)z^ + 1. On comparison with (I5])-(|6]), we arrive at the 
desired relationship between A and ju. □ 

3 Formulas for calculating a gear mapping 

Our numerical work will rest heavily on the following fact. 

Proposition 3.1 Let D be a pregear symmetric in M. Suppose that D is not 
a gear. Let p-i,pi G dD (p_i < pi) be the extreme points of the real interval 
UnM (i.e., the midpoints of the non-tooth edges of D). Let p be the common 
radius of the tooth edges, and let p G dD be an interior point of either tooth 
edge. Let v be the inward pointing unit normal to the tooth edge at p. Define 


c = p + rv, d = — (Imc)^)^^^ , 6 ^=Rec±d, 



Figure 2: Structure of pregears as in the proof of Proposition 13.II Note that 
precisely one of is interior to the pregear. 





and 


z — b 
z -b+ 


ifp-i <b < pi, 


Th±{z) = < 

Then the image G = Tb±{D) 
parameters are given by 


- 6 + 


ifp-i < 6+ < pi. 


z — b~ 

is a gear with gear center at the origin. Its gear 


^ Tb±(pi) 

^ T,±(p_i)’ 

7 = argTfe±(p). 


Proof. By construction, c is the center of a tooth edge. By Proposition 12.11 
the tooth edges must meet, which is impossible unless | Imc| < p (Figure [2]). 
Thus d is real (we assume positive), and the intersection points of the tooth 
edges are thus b^, with b~ < b'^. The dehnition of Tb±{z) is arranged so that 
Tb± leaves MU{cx3} invariant, it has its pole exterior to D, and T^^{x) > 0 for 
a; G M. The points are sent by Tb± to 0, oo, and hence the images of the 
edges are straight, so G is a gear by Proposition 12.11 which further implies 
that the images of the remaining edges are circles centered at the origin of 
radii —Tb±{p-i) and Tb±{pi). Since Tb±{p) lies on the upper tooth edge and 
this edge prolongs to pass through the origin, the formulas stated for /3, 7 
hold. □ 

An operation equivalent to changing the 2-jet at the origin may be effected 
with the aid of a self-mapping Tq of D dehned by (jS]). We use the following, 
which is verihed by a direct application of the Chain Rule. 


Proposition 3.2 Let f and g be two mappings such that Jf{zo) = (oq, oi, 02 ) 
and Jg{ao) = (&o,^i,^ 2 )- Then 

Jgof{zo) = {bQ,aibi, a ^62 + fl2&i)- 


In particular, 

JforM = (1-g^) ((l-g^)/"(-g)+ 2g/'(-g))) . 

The mapping f o Tg does not in general £x the origin. To achieve this 
condition we must use h = Tb± o f oTg where Tb± was defined in Proposition 
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EH If p_i < b < Pi, then the 2-jet of the Mobius transformation Th± at b 

is 


^ (°’ 6 +- 6 -’ {b+-b-y') ■ 

Otherwise, if p_i < b'^ < pi, then the 2-jet of the Mobius transformation 
at 6+ is 


^T.n6+) = (o, 


6 + — 6 “ ’ ( 6 + — b~Y 
We now apply Proposition 13.21 to obtain 


A(0) = 0, 




b+ — b 

_^-q 

b+ 


/'(-?), 


±2(1 - g2) 


2\2 




(6+ — b~Y 

^((1 - q‘)f"(.-q) + 2 «/'(-«)) 


(16) 


and as a result we have the following full description of the normalized gear 
mapping. 


Proposition 3.3 Given t,X, let f denote the solution of Sf = Rt^\ nor¬ 
malized by J/(0) = (0,1,0). Suppose that the image /(D) is a pregear. Let 
q = —f~^{b~) or q = —f~^{b'^) where the pregear center b~ E M. or b^ E M. 
is as described in Proposition \3.1{ Then the solution h of the Schwarzian 
differential equation Sh = {Rtp ° Tq){T'^‘^ normalized by the 2-jet of l[T5\) is 
a gear (not only a pregear) mapping satisfying h{0) = 0, h'(0) > 0. 

In order to apply Proposition 13.31 it is necessary to know the radius p in 
Proposition 13.11 The following result will hll this need. 


Proposition 3.4 Let f be holomorphic near zq = f'izo) 7^ 0, and sup¬ 
pose that for |;2| = 1 near f{z) lies in an arc of a circle of some radius 
p. Then the curvature of this arc is 


1 

P 


\fU)\ 


^1 ± ^0 


fM] 

fi^o) J ' 


A simpler version of this result was used in [5l |T6] . The choice of sign in this 
curvature formula corresponds to whether f{rzQ) enters or leaves the circle 
containing the arc when r < 1 increases to r > 1. 
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Proof. It is sufficient to consider the special case that Zq = 1 (and then 
apply the result to /(e®‘°z)). Let c denote the center of the circle \w — c\ = p 
containing the image arc. By the Reflection Principle, 


(fix) -c){f{l/x) -c) = 
for 1 — e < X < 1 + e. From the derivative 

- c) - \{f{x) - c)f'{l/x) = 0 

Ob 

we see that /'(1)(/(1) — c) = (/(I) — c)/'(l), or equivalently 

Im/'(1)(/(1) - c) = 0, arg(/(l)-c) = arg/'(l) modvr. (16) 
Take a further derivative and evaluate at x = 1: 

/"(i)(7(Ty^) - 2ir(i)p+ 2 f(i)(/(i) - c)+m(/(i) - c) = o, 

from which it follows that 

2Rem(/(1) - c) - + 2Re7^(/(1) - c) = 0 


Re((r(l) + /'(l))(/(l)-c)) = |/'(l)r. 

As a consequence of flT^ . 

Re(/"(l) + /'(l))(/(l)-c) = Re (^(/"(l) + /'(l))|/(l)-c|L(^) 

= ±|/(l)-e|Re(^(/"(l) + /'(l))lLMj, 

We therefore have 

±1 Re ((/"(1) + /'(1))^) 1 p A I 

1/(1)-ci if(i)p If (1)1 " V /'(!);■ 
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Numerical computation of gear parameters 


4.1 Auxiliary second-order linear ODE 


The basic facts relating the second-order linear differential equation 2y" + 
Rt,\y = 0 to the Schwarzian derivative of a conformal mapping are widely 
known [20] and have been applied in many contexts. A conformal / mapping 
with Schwarzian derivative Rt^\ is obtained as a quotient 

/ = - (17) 

yi 

of any two linearly independent solutions of the differential equation 


2y" + Rt,xy = 0 (18) 

or equivalently, as an antiderivative of yi^ (we assume yi nonvanishing). 

Thus given parameters t, A, it is a straightforward matter to compute / 
numerically such that 

Rf = Rt,\- (19) 

A particular solution will depend on the normalization chosen for yi, y 2 at a 
given base point. In general the image D = /(D) will not be a gear. Assuming 
it is at least a pregear, in order to apply Proposition 13.41 and Proposition 13.1[ 
we must be prepared to calculate the 2-jet J/ at certain boundary points, 
which may be accomplished by the relations 


= A f 

yi 




where we assume that the constant yiy '2 — y 2 y'i is equal to 1. 


( 20 ) 


4.2 Computational procedures for the disk 

We describe here some aspects of the calculation for mappings defined in the 
disk D. Similar considerations apply to mappings defined in the rectangle 
cis well as further observations we will make in 14.31 First we discuss 
calculation along radii. 
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Consider a fixed value of to- When we write equation 0181) by parametriz¬ 
ing along the radius from 0 to Zq = 6**“, it takes the form 


r]"{r) + Rt^x{re''^°)ri{r) = 0 


( 21 ) 


with ri{r) = |/(re**°). We use to = 0,7r/2,7r, i.e., we calculate the values of / 
along the rays from 0 to 1, t, —1. We consider two general approaches to the 
computation. 

4.2.1 Calculation via the 2-jet. 

The hrst, direct approach is simply to solve the equations numerically 
along radii. The values of the solutions at r = 1 together with relation (|20D 
and Proposition 13.41 and Proposition 13.11 permit us to calculate the transfor¬ 
mation T{z) which sends the pregear to a gear. Then F = T o / is a gear 
mapping with Sf = Rt^\. 

4.2.2 Calculation via the spectral parameter power series. 

In the second approach, using ([5])-([6]) we rewrite flT8|) as 


y" + i’o.ty = 


( 22 ) 


and again integrate along the three rays from 0 to 1, i, —1. We are following 
very closely the notation used in [5l [16], but now the functions 'ipo,'ipi are 
somewhat different since we are considering a different class of quadrilaterals. 
Then we make use of the SPPS representation of the solutions of fl22p . as was 
done in [S1|TB]. A summary of the main result establishing this representation 
is given in the Appendix. The procedure for gear mappings is as follows. 

• Calculate the SPPS integrals routinely in terms of the data 

"00,0 These are functions on 0 < r < 1 which depend on t but 

not on A. These indehnite integrals produce the coefficients of power 
series fl3^ in A dehning linearly independent solutions rii,rj 2 of fl22|) in 
D, normalized by 1-jets (?7i(0),77^(0)) and (772(0),772(0)) equal to (1,0) 
and (0,e**°), respectively. This normalization corresponds via fll7p to 
the conformal mapping / normalized by the 2-jet 


J;(0) = (0,1,0) 
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at z = 0. 


• Evaluate the coefficients of the power series at r = 1. This 

defines ?72(1) as functions of A (separately for each of the three 

rays). 

• Express the 2-jet Jf as a function of A by means of fl2U]) . 

Some variants of these approaches will be mentioned later. 


4.3 Computational procedure for the rectangle 

4.3.1 Condition to make pregear into a gear 


Analogously to the mappings of the disk, the Sg can be expressed in terms 
of solutions of the ordinary differential equation 

2/(0 + ^rAOviO = 0. (23) 

Two particular solutions are normalized by Jy^^O) = (1,0), Jy2{0) = (0,1). 
The quotient go = 1/2/1/1, normalized by Oo(0) = (0,1, 0), generally does not 
have a gear as image; the gear mapping must be sought within the general 
solution of Sg = i-e. 

^ ( 24 ) 

" cf^+d dyi + cu,_ ' ' 

with ad — be = 1. Thus 


9(0) = 


b 

c 




ad — be 
6/2 


/(O) 


—2 {ad — be)e 


We may assume a,bicA ^ because we are only interested in the case that 
g'{C) G M for C G M. We may further suppose that (7(0) = 0, (7'(0)=1 without 
affecting whether the image is a gear or not, so with this simplification we 
have 6 = 0 and a = d. Therefore we write 


9 = 


92 

91 + C192 


(25) 


where a = c/a G M, and the question can be restated as how to choose a to 
make the pregear a gear. 
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Proposition 4.1 Let the mapping g be determined by where yi, y 2 are 
normalized solutions of / d51) in Rq. Suppose that the image (7o(-Ro) of go = 
112 /yi is a pregear. Then g is a gear mapping if and only if a is a root of the 
guadratic eguation 

(lmj/21/2 U3/2) + (Im + y 2 y[) 1^3/2) a + lmyiy[ 1^3/2 = 0 

and the image g(JD)) is bounded. If the guadratic eguation has no roots, then 
go does not even produce a pregear. 


Proof. By Proposition 12.11 the condition for g{Ro) to be a gear is that the 
image of the upper edge of dRo should be straight. Parametrize this edge by 
0 J 3/2 — s for s G [0, cji). Since generally yi, y 2 do not take real values here, we 
introduce an auxiliary function n(s) = e~^'^y{u 3/2 — s) as the general solution 
to 


2v" + ip {—— s)v = 0. 


This differential equation has real coefficients, and the initial values are 


0(0) = “'(o) = 


where y is the general solution to fl2^ . If we can find 7 so that both n(0) 
and n'(0) are real, then we have v{t) real for all t. Clearly the condition is 
that y'{uJo/2) is a real multiple of yiuJo/2). Further, when this holds, we have 
y{uJo/2 — s) = e^'^v{s), so aigy{uJo/2 — s) is constant. Therefore argl/y^ is 
constant on the top edge of dRo- 

From fl2^ we have that g' = (j/i + aj/2)“^. Writing y = yi + ay 2 , our 
problem is to choose a so that 


y[ + ay 2 is a real multiple of yi + ay 2 at ( 


2 ■ 


The statement of the Proposition now follows from the fact that given any 
zi, Z 2 , wi, tC2 £ C, a G M, then wi + aw 2 is a real multiple of zi + az 2 if and 
only if 

(Im Z 2 W 2 )a‘^ + (Im {ziW 2 + Z 2 Wi))(^ + ZiWi = 0, 

as can be seen by multiplying the numerator of (wi + aw 2 )/{zi + az 2 ) by the 
conjugate of the denominator and then taking the imaginary part. □ 
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4.3.2 Computational procedure 


Now we can specify the numerical implementation of Proposition 14.11 While 
(123]) holds for ( & Rq, numerically we solve it first along the real interval 
[ 0 ,a;i/ 2 ] and evaluate the 1 -jets of yi, y 2 the endpoint to obtain 

(y) = (h.K), Jm (y) = ('’2' '’j)' (26) 

Then we solve the initial value problem 2u'\s) — y){uJi/2 + is)u{s) = 0, 
'^•ui(O) = ( 1 , 0 ), ^^ 2 ( 0 ) = ( 0 , 1 ) for t G [ 0 , |a; 2 |/ 2 ], obtaining at the endpoint 

=("2><=2)- (27) 

It follows from this that yi{ui/2 + is) = biUi{s) + ib[u 2 {s), y2{uJi/2 + is) = 
b 2 Ui{s) + ^ 62 ^ 2 ( 5 ), because the left and right sides satisfy identical initial 
conditions. Evaluate this at s = \uj2\/2 using (123]) . (127]) to obtain 

4i = (& 1 C 1 + f6'iC2, fe'ic'a - i6ic(), 

= (&2C1+ f62C2, &2C2-*^ 2 c'i). ( 28 ) 

The values 61 , 62 , ^ 2 ; c(, C 2 , C 2 are all real. They give in (123]) the values 

required in Proposition 14.11 which can thus be found by two integrations on 

a real interval and two integrations on the right edge of dR^. 


5 Numerical results 

We give here some specific results of the numerical application of the opera¬ 
tions described in the previous section. 

5.1 Computing [3 and 7 as functions of t and A 

Radial integration. The hrst, most direct method discussed in 14.2.11 above is 
to substitute into ( 120 ]) the boundary values obtained by solving ( 12 T]) along 
three specihc radii in order to apply Proposition 13.41 and Proposition 13.11 
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Going back and solving (jl]) numerically on any radius, we calculate f{z) for 
any z G In particular, we have all boundary points of the image gear 
/(D) as long as we stay away from the prevertices (where Rt,\ has poles). 
Then it is easy to hnd the slope of the tooth edges, the gear center, and the 
extreme points /(—I), /(I) of the other two edges. From this we have a way 
of calculating /?, 7 as a function of t, A. This method of calculation was used 
to produce Figure [3], which suggested Conjecture 3.5 of [B] to the effect that 
for hxed /9, there are precisely two values of 7 for each t below a threshold 
value tp. 

SPPS integrals. The approach described in 14.2.21 was carried out in many 
examples. To apply Proposition lA.ll it is necessary to have a nonvanishing 
solution of the differential equation. This amounts to having a good initial 
guess for Aoo, which is a simple matter in the light of formulas fl29|) below. 

Once this procedure is carried out, we again obtain the data of Proposition 
13.41 and Proposition l3.ll but now expressed as a function of A. This permits to 
study the behavior of (3 and 7 as functions of A (for hxed t), and also to solve 
for A with desired properties. In particular. Figure H] gives an illustration of 
the curvature k(A) of the tooth edges for a hxed value of t. When k(A) = 0, 
the mapping / with Sf = Rt^\ sends the arcs [ 2 : 1 , 2 ^ 2 ] and [ 2 : 3 , 2 : 4 ] to subsets 
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Figure 3: Pairs ( 7 , /3) for hxed t = tt/ n with n = 3,4,..., 10. Larger dashing 
indicates larger value of t. 
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of straight lines. However, in general the solution / is not univalent: one 
must choose here the largest value of A. It is interesting that by this method, 
one can solve for any desired value (not necessarily zero) and thus obtain 
pregears with prescribed radii for the tooth edges. 

As described in 14.2.21 we may substitute the SPPS formulas into those of 
Proposition 13.11 to obtain the coefficients of the Mobius transformation T{z) 
as functions of A. Applying this “symbolic” T to the SPPS series for p_i(A), 
Pi (A), p(A), we obtain formulas /5(A) and 7 (A) as combinations of power 
series (one could apply the Cauchy rule for products of series, together with 
inversion of series, to obtain a single power series in each case, but this is is 
rather complicated and fortunately is not necessary). An example is shown 
in Figure [5l When combined, ( 7 (A),/5(A)) parametrize a graph such as in 
Figure [3l 

Inasmuch as this calculation depends on the values of Rt^x{z) for 2 ; on the 
radius from z = 0 to Zo = i it necessarily involves complex arithmetic. As an 
alternative to this method, we can also use the other criterion of Proposition 
O that the non-tooth edges be concentric. By integrations along the real 
axis, which only involve real values of Rt,\, we can calculate the centers of the 



Figure 4: Graphs of the curvature k{X) of the tooth edges of the family of 
pregears corresponding to Rt\ for t/n = 0.1, 0.2, 0,3, 0,4. 
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Figure 5: Graphs of /9(A) and 7 (A) produced by SPPS formulas, for f = vr/d. 

circles containing these edges as /(—I) + 1 /k_i and /(I) + I/ki, where k_i 
and Ki are the corresponding curvatures. Subtracting the SPPS formulas for 
the values of the centers produces a function of A which vanishes when the 
centers coincide. Numerical experiments indicate that either approach seems 
to work equally well. 

Repositioning of gear center. To avoid possible confusion we reiterate that 
with the condition J/(0) = (0,1,0), the value /(O) is not likely to be the 
gear center wq of /(D). Therefore such mappings are quite different from 
the solutions F of the classical integral representation of gear mappings as 
in p[iniE 2 ]. Consider a general conhguration of prevertices 
(recall the discussion at the beginning of I2.ip . and let / be the solution of 
Sf = normalized by J/(0) = (0,1,0). Having calculated the data 

of Proposition 13.11 we can hnd wo, and presumably having already solved 
/ = y 2 /yi along [— 1 , 1 ], it is possible now to approximate p = f~^{wo) 
by numerical inversion. Then the composition F = f o T^p (which has a 
different 2-jet at the origin) sends 0 to Wq. Further, F(D) = /(D). Clearly 
Sp = Rt',t' ,x where T_p(e*h) = T_p(e*h) = Based on numerical 
examples of this procedure we are led to conjecture the following. 

Conjecture 5.1 Given ti,t 2 , 0 < H < ^2 < 7 r/ 2 , there is a unigue A G M 
such that the solution f of Sf = normalized by J/( 0 ) = ( 0 , 1 , 0 ) is a 

gear mapping. 

This is equivalent to the statement that if there is a conformal mapping of 
gears —)■ respecting the vertices and with 2 -jet of the form ( 0 , r, 0 ) 

at the origin (r > 0), then (3 = (3' and 7 = 7 '. 
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Integration in a rectangle. The method described in 14.3.21 is easily applied, 
integrating from 0 to Ci;i /2 and from Ui/2 to (wi + u; 2)/2 to hnd the required 
parameter a for converting the gear to a pregear. In the hrst integration 
one may save the values of yi,y 2 at points along [ 0 ,a;i] and then use them 
as initial values for integrating upwards or downward on vertical segments 
passing through those points. An example is shown in Figured The image at 
the right results from the value of a which produces an unbounded domain. 
The complement of this unbounded domain is again a gear domain; we know 
of no relation between this and the bounded image. 




Figure 6 : Mapping of rectangle to gear (r = 1.5). 


5.2 Inverse problem: prescribed gear parameters 

As occurs in many contexts in conformal mapping, the interesting problem 
is to hnd the auxiliary parameters which produce a given geometry. First 
we consider the following simpler question: Given t, determine A so that the 
image of the gear mapping ft x has gear ratio /3 (or alternatively, gear angle 
7 ). Since t is hxed, we can calculate the formula approximating k(A), solve 
k{X) = 0, and then having t, A we determine (3. This can be repeated as 
necessary, by a process of successive approximations, to make (3 have the 
desired value within specihed accuracy. We have carried out this approach 
successfully. A more direct method is to use the SPPS formulas for /5(A), 
7 (A) (recall Figure E]), and solve them directly for the desired values of / 5 , 7 . 
However, for extremely small values of t this does not give good results unless 
a great number of powers are taken in the SPPS formulas. 

The method of Broyden [S] may also be used to hnd zeros of the mapping 
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(t, A) ^ (/9 — (3o, 7 — 7 o), without recourse to the SPPS formulas. This turns 
out to be extremely fast (thousands of solutions in less than a second on 
an ordinary portable computer). However, the Broyden method works by 
jumping around unpredictably in from the initial guess for (t, A), and 
may fail by leaving the (t, A) region where (/ 3 , 7 ) is well dehned. Figure 
[7] shows values where this method succeeds starting from the initial guess 
(t, A) = ( 7 r/ 4 , 0 ). We will not pursue further the question of improving the 
initial guess. 

6 Internal structure of the region of gearlike- 
ness 

In [ 6 ] it was shown that the region of gearlikeness in the {t, A)-plane is 


y 

3.0 - 

2.5 - 

2.0 - 

1.5 - 

1.0 - 

0.5 - 

0.0 - 

0 2 4 6 8 10 

Figure 7: Experimental region of success of Broyden’s method for inverting 
the parameter correspondence, applied for a rectangular grid of 1.1 < /5 < 
10.0 and 0.1 < 7 < tt — 0.1. 


6 = {((,A): A,-<A<A+}, 
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where 


A 



We use our numerical methods to obtain a very illuminating picture of the 
structure of Q as related to {(3,'j). Observe that the vertical cross-sections 
of Q (representing gears of a given conformal modulus M{t)) are of common 
height = 1/2 for all t, and 



An = lim A 


We use the methods of the previous section for {(3, 7 ) 1-4 (t, A) to calculate 
the level curves of the geometric parameters. Figure [S] shows the subsets of 
Q of constant log/3 = 0.2, 0.4, 0.6, 0.8, 1.0; 1.25, 1.5, 2.0... and of constant 
7 = O.lvr, 0.27r, ..., O.Ovr. 

Individually, the /3-curves accumulate only at the two extreme points 
(0, Aq ) and (0,Ao), never at interior points (t, Af) of the lower and upper 
boundaries. At these extreme points /3 —)■ 1 while 7 —)■ vr, 0 respectively, and 
the gear degenerates to a disk in either case. It may appear paradoxical 
that the limiting Schwarzian derivatives {z) are not identically zero; 
however, the pullbacks of the Schwarzian derivatives according to Proposition 
13.31 by appropriate Tq (with q depending on A) do vanish in the limit. For 
hxed 7 , as f —)■ 0 or t —)■ 7 r /2 we have /3 —)■ cx) or /3 1 respectively. For 
hxed /3, as t —)■ 0 we have A —)■ Ag or A —)■ Aq and then 7 —)■ vr or 7 —)■ 0. 

This is one of the qualitative results on the conformal module proved in [ 6 ]. 

Each level curve for fixed 7 intersects the A-axis in a value, which might 
be termed limt^g corresponding to a degenerate gear. As /3 —)■ 00 for 
hxed 7 , the gear tends to the union of D with a full sector of angle 27 . 
Letting ti —)■ 0 in the gear mapping, we are led to apply the general formula 
(4) to this circular triangle with angles 37 r/ 2 , 37 r/ 2 , 7 and prevertices 
1 , to obtain 



(27/7r)^(cost2 - 1) 


(z — l)^(z^ — (2 cos 3 : 2 ) 2 : + 1) 

(cost 2 — 1)((52:^ — 14^ -|- 5) cost 2 + 7z^ — lOz + 7) 


2{z — iy{z‘^ — (2 cost 2 )z -I- ly 
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Figure 8 : Region of gearlikeness Q foliated by (3 level curves (left) and 7 level 
curves (right). 


whereas 


8A(cost2 — 1) 

{z — — (2 COSt2)z + 1 ) 


5 sin^ t2 

2{z — iy{z‘^ — (2 cost 2 )z + 1)2 


Equating these two Schwarzian derivatives we hud that 


lim 

t —>^0 



(30) 


independently of ^ 2 - This formula is conhrmed by Figure [S] (right). 

Goodman noted in [10] that when one limits the discussion to a particular 
7 (and Exes the normalization as /(O) = 0 = f'{0) — 1), a relation is deter¬ 
mined between the parameters ti and ^ 2 , stating that he could not calculate 
it except for the particular case 7 = 7r/2, where he found that 

1 

costi = 1 , COSt 2 = 2 ’ 

i.e. = 0 , /5 = cxo as the outer vertices of the gear have coalesced at 00 , 
while ^2 = 7 r/ 3 . Using Goodman’s explicit mapping formula 

4 2(1-z +z'^)-2 + 3z + 3z^-2z^ 

~ ^ z(l-z) 

and comparing Sf with Ro,tt/ 3 ,x, one Ends readily that A = 0, which is thus 
the value of limi_j.o conErming fl5U]l for this case. 
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7 Applications and conclusions 

We close with some brief applications of our results on gear mappings. 


7.1 The first Maclaurin series coefficient 

In [10] Goodman left unsolved the problem of calculating the ratio 6 i//(l), 
where hi = f'{0) is the hrst Maclaurin coefficient of the conformal mapping 
h of © onto a one-tooth gear domain that maps the origin to the gear center. 
This ratio was expressed in terms of singular integrals in [22] and later in [^, 
where it is shown that 



The integrand presents rather complicated singularities at the endoints of 
integration. 

In the construction of Proposition lTSl we obtained the 2-jet J/(0) via fflSj) . 
so we have f'{0). It is also a simple matter to calculate /(I) numerically. 
Thus one may evaluate integrals of the form given above via solutions of gear 
mapping problems. 

7.2 Module of the complement of annular rectangle 

Consider the bounded region A = with boundary dA = {e*^: 'y < 9 < 
27r — 7 }U{re*'^: 1 < r < 7 < 0 < 27r — 7 }U{re“*'>': 1 < r < 

i.e., an “annular rectangle” as in Figure |9] subtending an angle 2 ( 7 r— 7 ) within 
an annulus of radii of ratio f3‘^. While it is easy to map a disk or rectangle 
to A, it is not known (cf. [T71 p. 122 ]) how to obtain an expression for the 
conformal module of the exterior A* of A in closed form in terms of of the 
modulus of A. (This is the situation with quadrilaterals in general. There has 
been a surge of interest recently [H [HI [ 121121 ] in the question of numerical 
calculation of exterior moduli of topological quadrilaterals.) Here we give a 
numerical solution to this problem. 

Theorem 7.1 Let A be the annular rectangle with inner radius 1 and outer 
radius and angle 2{n — 7). The conformal module of the exterior A* = 
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Figure 9: Complement of annular rectangle is formed of a gear and its reflec¬ 
tion along its B-arc. 

(C U {c)o}) \clA is half the module of the gear domain with gear angle 7 and 
gear ratio {3. 

Proof. Consider the conformal mapping g: Ro —)■ of Proposition 12.21 

The Schwarz reflection applied across the right vertical edge of Rq produces 
a mapping onto A* from a rectangle having double the width of i?o- hi 

Given an annular rectangle the values t and A corresponding to /3 
and 7 can be computed as described in Section I5.21 This means that ft^x 
is a mapping onto a gear with conformal module M{t). According to 
Theorem 17.11 the module of Af^^.y is therefore calculated numerically to be 
M{t)/2 . 


7.3 Multitooth gears 


Let /: D —)■ be a conformal mapping to a normalized gear domain, 

/(O) = 0. The function fn{z) = x/fi^ is a mapping to a regular n-toothed 
gear. Applying the Chain Rule to o /„ = / o where Pn{z) = z"', and 
using iSp„ = (1 — n‘^)z~‘^/2 we find 


i-n^gz? 

2 fn{zf 




n 


Since both Sf and are rational functions, it follows that {fUfn)'^ is a 
rational function. With somewhat more work one recovers the formula of 
[To] of which equation (2) of [6] is a particular case. 

The image fn{Gg^.f) has gear ratio and each tooth subtends an angle 
oi2'^/n with spacing of 2 ( 7 r — 7 )/^ between consecutive teeth. By means of 
these facts it is simple to calculate any desired regular multitooth domain. 
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Figure 10: Single-tooth gear (left) calculated to generate a multitooth version 
of prescribed geometry (center). Detail (right) includes level curves corre¬ 
sponding to \z\ = 0.9, 0.91,..., 0.99, 0.991,..., .999, 0.9991,..., 0.9999. 

We illustrate with the example of n = 10 teeth, having intertooth space equal 
to the tooth width, and gear ratio arbitrarily chosen as 1.3. Thus we have 
(3 = 1.3^° ~ 13.79, 7 ~ 7r/2. By the methods given in Section [521 we hnd 
that the one-tooth gear with these parameters is obtained by f ~ 0.6024, 
A ~ —0.0029. The result is in Figure [TOl We stress that the renormalization 
worked out in Proposition 13.3! is essential, since without the condition that 
/(O) = 0 is the gear center, the application of Pn will not work. 

It may be noted that when the number n of teeth is large, the circu¬ 
lar arcs are approximated by straight lines, and the gear mapping may be 
approximated by a Schwarz=Christoffel integral. 


A Appendix: SPPS method 

The sequence of iterated integrals generated by an arbitrary pair of func¬ 
tions (go, gi) is dehned recursively by setting Jq = 1 identically and forn > 1, 

In{z)= r In-l{C)qn-l{OdC (31) 

Jo 

where g„+ 2 j = g^ for j = 1, 2,... 
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Proposition A.l [TS] Let 'ipQ and 'ipi be given, and suppose that is a 
nonvanishing solution of 


y'L + ^0 2/oo = AooV'i Voo 

on the interval [0,1], where Aoo is any constant. Choose go = = 

fjiylo and define to be the two seguences of iterated integrals 

generated by (go,<?i) by {qi,qo), respectively. Then for each A G C the 
functions 


yi = yooJ2{X-Xoo)^x(^^\ 

k=0 
oo 

1/2 = yooJ2{X-Xo.fx^^^+^'> (32) 

k=0 

are linearly independent solutions of the equation 

y" + fioy = Xfiiy (33) 

on [0,1]. Further, the series for yi and j /2 converge uniformly on [0,1] for 
every A. 

It is a straightforward matter to obtain the appropriate linear combina¬ 
tion of solntions with desired 2-jets at 2 ; = 0, for example (0,1) and (1, 0). 
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